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CN ■ Abstract 

j'Qi In applied mathematics generally and fluid dynamics in particular, the role of complex variable methods 

' is normally confined to two-dimensional motion and the association of points with complex numbers via the 

assignment w = x + iy. In this framework 2D potential flow can be treated through the use of holomorphic 
functions and biharmonic flow through a simple, but superficially non-holomorphic extension. This paper 
explains how to elevate the use of complex methods to three dimensions, using Penrose's theory of twistors 
as adapted to intrinsically 3D and non-relativistic problems by Hitchin. We first summarize the equations 
of 3D steady viscous fiuid fiow in their basic geometric form. We then explain the theory of twistors for 3D, 
S I resulting in complex holomorphic representations of solutions to harmonic and biharmonic problems. It is 

^> - shown how this intrinsically holomorphic 3D approach reduces naturally to the well-known 2D situations 

•^J ' when there is translational or rotational symmetry, and an example is given. We also show how the case 

1—^ , of small but finite Reynolds number can be integrated by complex variable techniques in two dimensions, 

j^ • albeit under strong assumptions. 

C/3 ' Key Words: Fluid dynamics, Twistor theory. Three dimensions, Stokes flow. Complex variables. Stream 

function, Biharmonic equation. Low Reynolds number 
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>>; 
^ '. 1 Introduction 

The use of complex variable techniques in applied mathematics, and especially fluid dynamics, is dominated 

by two-dimensional applications through the prescription 
> , 
\l , w — X + ill . (1) 

On the other hand, in the context of relativistic physics in four or more dimensions, the use of twistor methods 
due to R. Penrose and co-workers is becoming an ever more present tool in the hands of theoretical physicists. 
Ji]^ I Many of the key earlier developments are summarized by Penrose & Rindler (1984). A collection of more 

^"^ ■ recent developments is available (Mason & Hughston, 1990, Mason et al, 1990) and some state of the art 

applications to superstring theory are given in a very recent thread of work initiated by Witten (2003) (see 



o 



also the thread on the web site |http : //arxiv . org/cits/hep-th/0312171) 



. , ^ , In the important intermediate case of non-relativistic problems in three spatial dimensions, rather than 

j^ ' working with a copy of C for the two-dimensional case, or with CP^ (the twistor space for Minkowski space- 

JH I time) we work with the twistor space for three dimensions: TCP^ - the tangent space for the Riemann 

sphere. This case has not received quite such wide attention. A notable exception is the paper by N. Hitchin 
(1982), where a summary of the basic theory appeared as a preparation for an application to monopoles. In 
a notable Appendix, Hitchin (1982) also explained the use of TCP^ to explain the geometrical underpinnings 
of Weierstrass's holomorphic solution of Plateau's problem (minimal surfaces, or soap bubbles) (Weierstrass, 
1866). Hitchin lectured in Oxford and Stony Brook on the methods in the 1980s, but those talks have not been 
published. A summary of the some of basic ideas accessible to non-specialists is now available (Shaw, 2006), 
based on developing the three-dimensional theory via the relativistic method, and including applications to 
Laplace's equation in three dimensions. These methods give the complex geometrical underpinnings of contour 
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integral formulae obtained by more ad hoc methods and known for some time - see for example Whittaker's 
work (Whittaker, 1903). The applications of this approach to various types of problem with axis-symmetry 
has been given by Mason (1990) and Ward (1981). 

The focus of much of the published work has been on time- independent problems within the general context 
of theoretical relativistic physics. In this paper the idea is to present such methods as being a routinely useful 
tool in traditional applied mathematics. To this end, an example of the application of twistor theory to 
viscous fluid flow is presented. In particular, the solution of various biharmonic problems will be presented 
using contour integral techniques. The ultimate goal of this work is a better understanding of the Navier- 
Stokes equations through the geometry of holomorphic complex variable techniques. At first sight, even our 
most basic goal might seem to be an unreasonable proposal. For example, the biharmonic equation in two 
dimensions, with the w ~ x + iy prescription, amounts to 

5^5S* = (2) 

with the general real solution 

^ = ^wf{w)+giw)} (3) 

where / and g are both locally holomorphic. This is generally regarded as going outside the holomorphic 
context as it involves w in an essential way. We shall show that equation (1.3) is in fact the two-dimensional 
projection of an essentially holomorphic three-dimensional result. 

One goal of this work is a better understanding of the Navier-Stokes equations through the geometry of 
complex variable techniques. There is already progress in this area. A characterization of the full steady 
problem in 2D has been given by Legendre (1949) and Ranger (1991, 1994). Work on the Kahler Geometry 
associated with the Navier-Stokes problem has been given by Roulstonc et al (2005), as well as related char- 
acterizations Roulstone et al 2009. It would give considerable insight to draw together theses various threads. 
We cannot expect to be able to analyze the full Navier-Stokes system in three dimensions using twistor meth- 
ods if we cannot first treat the biharmonic limit, and this may well give insight into the full problem - this 
current paper suggests it is at least straightforward to treat the biharmonic case. 

There are also several ways of attacking the problem based on different choices of physical equation that 
we wish to transform into a holomorphic geometrical description. For example, given that the momentum 
equation is fundamental, one could work with the equation directly, or its curl, or its divergence. Here we work 
with the curl in the guise of the vorticity equation. Roulstone et al (2005) work with the divergence. Similarly, 
here we work with a vector potential for the velocity. One might also consider working with the Clebsch 
potential representation, as discussed, for example, by Milne-Thomson (see Section 21.021 of Milne-Thomson, 
1996). Generalizations of this approach are much favoured for the treatment of fluids by action principles (see, 
e.g. Morrison, 1998). We will take the vector potential approach rather than the "Clebsch" route here for 
three reasons. First, in the latter approach the vorticity is already a non-linear function of the basic variables. 
Second, the relationship between action principles and twistor methods is an uneasy one, since the power of 
twistor methods lies mainly (so far) in an understanding of how the differential equations governing a system 
may be solved exactly by the introduction of some holomorphic structure. Dynamical configurations that are 
"off-shell" (in physics terminology), i.e. are not solutions of the equations of motion in space or space-time, 
do not so far appear to have a straightforward twistor description, though the recent work by Witten (2003) 
makes use of action functionals in twistor space to describe particle interactions. Finally, in taking the vector 
potential approach, we are able to make contact, for systems with certain symmetry, with the simple stream 
function approach. 

A key outcome of the development given here is that one can show show that there is a simple holomorphic 
framework for three dimensions, within which we can treat the vector biharmonic equation in holomorphic 
terms. We also show how to recover not only Eq. (1.3) for systems with translational symmetry, but also a 
contour integral representation of the Stokes stream function for the axis-symmetric case. These two symmetric 
systems emerge from a common 3D holomorphic framework. 

Although we are working with a biharmonic system motivated by fluid dynamics, there are of course 
applications to elasticity with regard to both the biharmonic equation in general (see Note 93 of Love, (1927) 
and Howell et al., (2008) for a modern treatment) and specifically complex variable methods (Muskhclishvili, 
1949, England, 1971). The methods discussed in this paper can equally be applied to the three-dimensional 
Lame or Navier equation for steady configurations of incompressible clastic media under a potential body 
force. 



The plan of the work is as follows. In Section 2 a brief summary of viscous fluid dynamics is given. In 
Section 3 a self-contained introduction to twistors for three dimensions is given together with the basic ideas 
of how to solve the scalar Laplace and biharmonic equations. In Section 4 this is extended to the vector 
biharmonic case and the relationship to planar and axis-symmetric stream functions is elucidated. In Section 
5 we take a first look at the non-linear problem arising from the case of small but finite Reynolds number 
in two dimensions. It is shown how this may be integrated (almost explicitly) using entirely holomorphic 
methods, under strong assumptions on the nature of the fiow. 

2 Steady viscous incompressible flow 

A large class of fluids can be characterized by their density, p, a scalar field not presumed to be constant, 
and their dynamic viscosity /i. The flow is characterized by a velocity vector fleld v, and an associated scalar 
pressure field p. Conservation of mass is expressed by the continuity equation 

^+^-ipv)=0 (4) 
and the conservation of momentum is expressed by the Navier-Stokes equation^j 

p{-^ + w-V w) = - Vp + ^V^u (5) 
If the fluid is incompressible in the sense that p is a constant in both time and space, we have the condition: 

V.w^^O (6) 
To analyze matters further, we introduce the vorticity vector 

a; = V X t; (7) 

In the following discussion we demand incompressibility but allow for non-zero vorticity. Using simple identities 
from vector calculus the Navier-Stokes equations may then be recast in the form 

p(^-l^xcj) + V(p+-pi;2) = -^V xw. (8) 

Taking the curl of this, we arrive at the vorticity equation 

^ -hv.Vw-w.Vi; = i/V^w (9) 

where the kinematic viscosity v = \il p. 

2.1 The 'stream vector potential' 

Since the velocity field is divergence-free, we may introduce a vector potential ^ such that 

u = V X * (10) 

and furthermore we may choose it so that it is divergence free: 

V.*.= (11) 

In theoretical physics, notably electromagnetic theory, this is known as setting a gauge condition. The tradition 
in fiuid dynamics is to mainly use the vector potential only when it can be reduced to a single function using 
some type of symmetry. The resulting object is a stream function. For example, planar 2D flow is obtained 
by setting (and note that this automatically satisfies satisfies the gauge condition) 

* = -^(x,y)e, (12) 



^In Sections 2 and 3 of this article, V^ acting on vectors should be understood as the ordinary Laplacian acting on Cartesian 
components. The effects of a non-constant basis are dealt with in Section 4. 



We will work with the full vector form. First of all we note that under the assumption that VP satisfies V .^ = 

cj = -V^* (13) 

and the vorticity equation becomes, denoting ^ by ': 

V^^ =-l ((V2^).V )V X £ - ((V X £).V )(V2£.) + V^*! (14) 

or indeed as 

74,T/ — J V7 s^ UXl N^ lT/^ sx V72lTf^ I Y72iTi 



V** = - <^ V X ((V X *) X V^*) + V^£ ^ (15) 

This latter representation of the Navier-Stokes equations is well-known in the 2D planar case (see for example 
Ockcndon et al. 2003) where it reduces, in the time- independent case, to the equation 

VV^i^^t^ (16) 

2.2 The biharmonic limit 

For problems where it is possible to identify a natural length scale L and a natural speed t/, it is normal 
practice to perform a non-dimensionalization of the variables and introduce the Reynolds number R = UL/v. 
Then the steady-state vorticity equation becomes (after rescaling the independent variables suitably): 

V% = i?|v X ((V x£) X V2£)l (17) 

The old historical approach to the limiting case when _R — >■ is to take the view that the non-linearities may 
be ignored (provided the non-linear term is well behaved) and the time-independent Navier-Stokes equations 
reduce to 

Vl. = (18) 

which is the biharmonic limit, also known as Stokes flow. We shall focus on the solution of this equation 
by complex variable methods. It is now well known (see for example. Chapter 8 of Van Dyke (1964)) that 
the neglect of the non-linear terms can lead to inconsistencies, as is evidenced by the lack of any solution for 
asymptotically uniform two-dimensional flow past a cylinder. However, in attempting to construct a twistor 
description of fluid flow we must be able to at least solve the biharmonic equation. It is to this that we now 
turn. 

3 Twistor solutions of the Laplace and biharmonic equation 

It is very well known that the Laplace equation can be solved in terms of holomorphic functions in two 
dimensions. Among devotees of twistor methods, it is also well known that this can be carried out in three 
dimensions. Our purpose is to extend this idea to the biharmonic case. 

3.1 Twistor space 

The twistor solution of the Laplace equation in three dimensions may be discussed at several levels. In 
presenting this material the simplest and hopefully most comprehensible route shall be taken, but it is necessary 
to be somewhat formal in order to state results correctly. This presentation follows that given by Hitchin 
(1982]j. The twistor space associated with R"^ is first, as a real space, the set of oriented straight lines in R^. 
Relative to some origin O, let r denote the position vector of the point on a given line nearest to O. Then r 
is orthogonal to the direction of the line, which we denote by u with u.u = 1. So the set of oriented straight 
lines is the set 

TS^ = \{l,u) e M^ X S-^ I r.u = ol (19) 



^This presentation also draws on Hitchin's unpublished lectures, particularly with regard to the solution of Laplace's equation 
and understanding the 2D limit 



where our notation emphasizes that this is the tangent bundle for the unit sphere. This set is also naturally 
the tangent bundle to a complex manifold, where we regard S^ as the Riemann sphere CP^. This complex 
tangent bundle, TCP^ is the twistor space of interest. 

The next question is how to define a point in ordinary space in terms of some structure on TCP^7 A point 
may be regarded as the intersection of all straight lines through it. This means that a point is necessarily some 
vector field in TCP^ that is defined globally. To see the implications of this we introduce two open sets that 
cover <CP^. We can take coordinates for the sphere as C on one patch (covering everything except infinity), 
and C = 1/C, on another patch, covering everything except C = 0. Over each of these respective patches we can 
define coordinates for the tangent bundle as (77, (.) and (ry, C), where the relevant vector fields are, respectively 

"re "s « 

Consider now a holomorphic vector field. On the (. patch it can be written as 

/o(C)|^ (21) 

for some /o, and on the C patch, it can be written as 

/i(C)|^ (22) 

for some /i . On the intersection of the two patches equality of the two representations gives us 

/i(r')(-c')^ = /o(c)^ (23) 

If we make a Taylor series expansion of both functions, /i(C) = S^o'^nC": '^e deduce that the coefficients 
ajj vanish if n > 2. That is, the global vector fields must be of the form, for example on the 77 patch: 

77(C) =a + bC + ce (24) 

so that such quadratics are the only holomorphic vector fields, and these correspond to points of C"', parametrized 
in some way by (a, b, c). 

Further analysis of this system allows the identification of real points in M.^, and the construction of a 
natural metric. The points are real if and only if 

c^-a AND b = b (25) 

The induced metric is proportional to the discriminant of the quadratic, and we shall normalize matters such 
that 

ds'^ = dx'^ + dy^ + dz^ = -db^ - dadc (26) 

The metric for quadratics corresponding to real points is therefore: 

ds^ == dx^ + dy'^ + dz^ = -db^ + dada (27) 

If we pick our coordinate system such that the real part of a is x, we see that we can take the imaginary part 
of a to be ±y and set b = ±2z. The convention is to set: 

vAC) = i^ + ^y) + 2^C -{x- iy)e (28) 

Note that here and elsewhere we will carefully avoid ewer writing z ~ x + iyl Note that although this particular 
Cartesian representation does not make manifest the action of rotations on M.^ it will turn out to be vary useful 
to deal with the cases of translational and axial symmetry. This will be discussed in Sub-sections 2(d) and 
2(c). 



3.2 Solving the scalar Laplace equation 

We consider a function /(ry, () defined on twister space. This can then be thought of as restricted to the special 
global sections of twistor space represented by rjr (C) , and the C-dependence integrated out by integration over 
a contour C. We set: 

</-(!)- //(%(C),CK (29) 

Jc 

Then (j> satisfies the scalar Laplace's equation. To see this observe that 

^^^^^#^M1-C^)'0U. (30) 

^^^|P^-^'a-c^)'0i,^. (31) 

a'/(%(0.0 ,.,,taV| ™ 

and that adding these three expressions with k — 2 gives zero identically for any choice of /. Note that many 
different choices of / will give rise to the same 4>. Such choices differ by the additions of functions that are 
holomorphic inside or outside of C, so that one must pursue a cohomological approach in order to state a 
formal isomorphism between structures on twistor space and solutions of the Laplace equation. 

3.3 Solving the scalar biharmonic equation 

The extension of the results for the scalar Laplace equation to the scalar biharmonic case may be pursued 
on several levels. Our argument does not require any advanced knowledge of twistor theory, though it would 
be interesting to characterize this analysis in terms of the appropriate cohomology structures. We shall pose 
the following question. How do we modify the integrand /(?7r(C)jC)j say to some holomorphic function g, to 
arrange that V^g = but V^g ^ 0? 

Let's consider trying to build g from / by multiplying by some prefactor h{r, C.), so that 

ff = MLC)/(%(C),C) (33) 

Now 

V^g = V^hf = fV^h + hV^f + 2V /i.V / = fV^h + 2Vh.Vf (34) 

where the last simplification arises as / satisfies the Laplace equation. If we furthermore choose h to be linear 
in r matters simplify further and we have 

V2g = 2V/i.V/ (35) 

Let us set, w.l.o.g. (other than excluding h being constant, which gives us harmonic solutions already under- 
stood), h — u{().r, so that V_h — u(C)- We also note that 

Y/ = |^Vry=|^(l-C^^(l + C^),2C). (36) 

Putting this all together, we arrive at 

V'g = V'hf = 2u{C).{l - C'Al + C'),2C)|^ = ^'^^K)^^)!^ (^^) 

We can now see that V^g = identically, while V^q does not vanish unless rjy^^QiC) = 0. To see what is 
happening, we can now make matters more explicit. We let m(C) = (""1(01 '"2(C)j ^3(0): then 

u-L = ui {C)x + W3 (C)y + Us {C)z (38) 

and 

%(c)(C) = ("i(C) + m2(C)) + 2u3(C)C - («i(C) - ^MO)e (39) 



In terms of these variables the proposed integral representation for solutions of the 3D scalar biharmonic 
equation is just 



*- / rfC 

Ic 



xui(C) + yu2{C)) + 2U3(C) 



/(%(C),C) 



or indeed as (this will prove useful presently when looking at axis-symmetry): with w — x + iy: 



*-y/^ 



w5_(C) + w.g+(C) + 2zM3(C) 



.nvr_{o,c) 



(40) 



(41) 



where g± (C) = ui (C) ± iu2 (C) ) • Note that the introduction of w is only a parametrization convenience - it has 
nothing to do with the underlying holomorphic structure of twistor space. 

3.4 The scalar biharmonic problem in 2D 

Suppose we want no z-dcpcndence. We set U3 = and w = x + iy, so that 



*- / dC 

Jc 
We can write this in the equivalent form 



wg^iO+wg+iC) 



finziOX) 



(42) 



^ = w dUi (ry. (C) , C) + ^ / dQh (%(C) , C) (43) 

JC Jc 

Now consider the second term. This is w(j){x, y, z), where (/> is a solution of Laplace's equation and is just 



{x,y,z) = / dCf2{r)r{0,0 
Jc 



(44) 



We want this not to depend on z either. But this looks awkward given that %(C) = {x + iy) + 2z( — {x — iy)C'^. 
It is not so awkward as it looks. What we want is translation invariance for d>. Note that 



ix,y,z + h/2)= / rfC/2(%(C) + K,C) 
Jc 



The equation we need is 
This does not require that 
Instead we need 



(f){x,y,z + h/2) = 4){x,y,z) 

f2iVriO+hCX)^.f2{Vr.iOX) 

/2 (^. (C) + hC,0^f2 (%(C) ,C)+9oivX,h)~giirj,Ch) 



(45) 

(46) 
(47) 
(48) 



where go is holomorphic on and inside C and gi is likewise outside. Cauchy's theorem then gives us the desired 
result. In order to give a clear calculation. Let's take C to be unit circle, or to be deformable to the unit 
circle. Now differentiate w.r.t h then set h — 0. We obtain, for some G,, 



c|| = Go(r7,C)-Gi(7?,C) 

We integrate this w.r.t. 77 and divide by (. We obtain, for some Hi, 

Ho{f^,C) H,{rj,0 



/2 



c 



c 



(49) 



(50) 



and recall that Hq must be holomorphic inside C and Hi holomorphic outside. Now we evaluate the integral 
of Eq. (3.26) using calculus of residues. The first term in Eq. (3.32) is easy, and we get 



27rii/o(77r(0),0) ^K{w) 



(51) 



for some functfon K{w), giving a contribution to (f) of wK{w). When we calculate the contribution of the 
second term of Eq. (3.32) to the integral of Eq. (3.26), we make the transformation C ^^ C ^nd obtain an 
integrand that is a function of 77 = (x — iy) — 2z( — (x + iy)C'^. Taking the residue at ^ = gives a function 
oi w = X — iy, also to be multiplied by w. 

The other two terms in Eq. (3.25) may be treated similarly. We end up with four terms contributing to '5: 

* = wK2{w) + wK2{w) + wKi{w) + wki{w) (52) 

When ^ is real we must have Eq. (1.3). So the fully holomorphic picture in three dimensions projects, via the 
calculus of residues, to a two-dimensional picture and generates the familiar yet superficially non-holomorphic 
two-dimensional representation of solutions to biharmonic (and Laplace) equations. In three dimensions our 
functions are contour integrals. 

3.5 The axis-symmetric scalar problem 

This problem is one of considerable interest, though it is important to realize that this is a different problem 
from the case of vectorial axis-symmetric flow - that issue and the link to the Stokes stream function will be 
discussed later. We go back to the representation 



*4//^ 



w5_(C) + w.g+(C) + 2zM3(C) 



/MCC) (53) 



with w = X + iy. We can regard this as three pieces, where we discard the factors of a half: 

^^=w [ dCg-{Of{Vr_iOX) = w^- (54) 

Jc 

*+ = ID / dCg+iOfiVriO, C) = ^^+ (55) 

Jc 

^3^z f dC"3(C)/(%(C), C) = z^Po (56) 

Jc 

In order to develop axis-symmetric solutions, we need to understand the action of the group of rotations about 
the z-axis. We need to bear in mind the formula 

ri^w + 2zC- wC,'^ (57) 

with w = X + iy. Under a rotation about the z-axis, z ^)- z and w — >■ exp{i(f>)w. This is compatible with the 
action (77, () -^ exp(i0)(77, (). In seeking axis-symmetric solutions for '0±,o we need to arrange that 

rfC.g-./^cxp(-j0)dCg-/ (58) 

dCg+.f -^ exp{i(j))dCg+f (59) 

dCusf ^ dCu3f (60) 

To treat all of these situations together, we consider the case where dC,h{rjX) -^ exp{im(f>)dCh{riX)- To this 
end we consider the contour of integration to be the unit circle and consider a basic set 

If r]" 

^"'"^2^y^'^c^T^ ^^^^ 

where for ipQ, m — 0, and for (j)±, m = ±1. So our task now is to calculate 

By multiplying these by the relevant factors of w,iu,z for m — —1,1,0 we get an interesting set of axis- 
symmetric biharmonic functions. The functions ipn.m- themselves are now contour integral solutions of 
Laplace's equation. This is of course of interest in itself. 



Note that a clue to what is going to emerge is that if n = 1, 2, 3 . . . , it is easily seen that the integral for 
ipn,m will vanish if there is a zero coefficient of C~^: ^-nd this will happen ii m > n or m < —n. So with n > 
we are only interested in \m\ < n (this is a hint to what special functions might emerge). 

To evaluate this set we let y = since ipn,m{r,9,(j)) = e^"^'^ip{r,9,0). Then, we have, in spherical polar 
coordinates, 

V-™,™ = ^ / dCr-'(2 cos(0) + (i - C) sin(f?))" (63) 

Parametrizing the integral as C = e'*, we obtain 

^n,m = ^?r- / dte'-^'icosiO) - sin(0) sin(i))" (64) 

Zw J 

On consulting Gradstheyn and Rhyzik (1980) equation 8.711.2, and performing some manipulations, we see 
that, discarding normalizations, if n 7^ — 1. 

rr"P™(cos(0)) n = 0,l,2,..., 
'" \p^Pr(cos(0)) n = -k-l,k = l,2,3,... ^ ^ 

When n = —1 matters are quite subtle as the integral branches depending on the sign of z! A full treatment 
of this is rather beyond the scope of this paper but we note that in this case, 

1 r C™ 

^-l..m = TT- / dC ^„\ _., (66) 

2m J w + 2zQ — wQ-^ 

The quadratic in the denominator has two roots C,± given by 

— z ± r 



C± = z — , r — J x^ + y^ + z^ , w — X — iy . (67) 

—w 

These roots are located, using standard spherical polar coordinates, at 

C+ = -e'*tan(^), C- = e'^^ cot(^) (68) 

and we can write 

^'-■■■-2^/-'< ((C-J)(C-C-)) « 

The details of the global evaluation of this are lengthy. We note here that when z > 0, |C+| < 1 and when 
TO > 0, the single residue inside the unit circle gives a value of 

w{C+-C-) 2r 

so in particular we obtain the Coulomb field in the region z > when to, = 0. The reader is invited to explore 
the other cases. 

4 The Vector Theory 

We now turn back to the treatment of the full vector potential theory. We need to elucidate the relationship 
between the vector potential approach and the more traditional use of stream functions. In the case of planar 
2D flow, there is almost nothing remaining to do in the biharmonic case. If we set the vector potential to be 

^^-^ix,y)e, (71) 

the divergence-free condition is satisfied identically and we merely need that 

Q2 g2 \2 



and this is all taken care of by the representation of ([5^. Note that this well-known intrinsically 2D and 
superficially non-holomorphic result is now clearly understood in 3D holomorphic terms. 



4.1 Axis-symmetric vector flow 

This is traditionally modelled (Milne- Thomson, 1996) in terms of the Stokes stream function ^^(r, 0). The 
components of the velocity field are given by 

«e = -—-T^^rr (73) 



r^ sin(6') 86 ' rsin(6') dr 

What this representation is really telling us, as is made clear in modern fluid theory, is that the vector potential 
^ for the flow is given by 

as is revealed, together with the fact that ^ is divergence-free, by elementary calculations with the curl and 
div operator expressed in a spherical basis. A further elementary calculation in vector calculus shows that for 
an axis-symmetric function f(r, 0) 

V X (V X — ^e^) = —^AE\f)e^ (75) 

rsm(t^) ^ rsm(t^) ^ 



where the operator E^ is given by 



2 d^f sin(6i) d / 1 a/' 



The biharmonic condition may be expressed as the scalar PDE 

E^^s = (77) 

However, this representation in some ways obscures the underlying simplicity of the problem. To see why, we 
need to work with the problem in the full vector form, and, perhaps surprisingly, recast it in a Cartesian basis. 
In this way we can use the contour integral technology already developed for the scalar biharmonic problem. 
We write the basis vector e± in the form 

— ue^ + xe„, —ve^ + xe„, 1 

" - -" - ' r5R{^«fe, + el} (78) 



~^ v/a;2 + 2/2 rsin(6l) rsin(6l)' 

where w = x -\- iy Q& before. We could equally well write this down in terms of w. Now recall that the full 
vector potential is given in terms of the Stokes stream function by the relation 

* = ^^^0 (79) 

rsm(t^) ^ 



and if Vl^s is real we can write the vector potential as 



= sin2(0)' 



^=^^^7^^^^^+^.^ (80) 



The components of this with respect to a Cartesian basis must satisfy the scalar biharmonic equation, or 
indeed, as a special case, the Laplace equation. We now appeal to equation (3.36), where we note that ^_ is 
just a harmonic function. It follows that we can write the parts of ^5 that are biharmonic but not harmonic 
in the form 

v&s = r2 sin2(0)g(r, d) = ww^^ J dC^f (^) (81) 

for some complex function /. Here g is harmonic and axis-symmetric and can therefore be written in terms 
of the tpn,o functions given in equation (3.44) or indeed in terms of normal Legendre functions and powers of 
r via a Laurent expansion of / in the form 

1 /" ^^ 

^s - r' sin\e)g{r, 9) = ww— dQ ^ ^(z« + 2zC - w^T (82) 



n— — 00 
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We argue that these relations are the natural axis-symmetric versions of (3.34). Of course, in general, we need 
to add in harmonic components, just as in the 2D planar case where we can add to ^ any pair fci(w) and 
k2{w) of holomorphic and anti-holomorphic functions. To treat this we look again at the representation in a 
Cartesian basis, this time in the form: 

* = . 'l„. (- sin(0)e^ + cos(0)e ) (83) 



We deduce that the function 



rsm{9) 



rsm(6') 

must be harmonic and therefore a solution of Laplace's equation with m = ±1 as described above. By 
packaging this up as before, we can write these harmonic contributions to ^5, say ^5^ in the elegant form 



l^/^<<?)-^/4<? 



•ls„ = ^ld(p{:i)+^l,K^y{:^) (85) 



for some choice of complex functions /3 and 7. 

We are finally able to give the proposed contour integral solution for the Stokes stream function for axis- 
symmetric biharmonic flow, by combining these expressions into the form: 



», , ^S,_l. / rf<l/| :i ) + ^ / <iC;5( ^ ) + ^ / d<A,( 1 ) (86) 



where 77 is written in terms of x, y, z and where /, /3, 7 have Laurent series expansions that generate expansions 
in terms of powers of r and regular (/) and modified (/?, 7) functions. 

4.2 A simple example to check it all works 

We can see that a choice of / constant gives a contribution to ^5 proportional to 

ww = r2sin2(6i) (87) 

We also know (at least locally) that the choice f{z) ~ 1/z gives a Coulomb field and a contribution to ^5 
proportional to 

^^rsm\e) (88) 

r 

Another interesting contribution can be generated by the choice I3{z) ~ 1/z^, where an elementary exercise in 
the calculus of residues leads to a contribution to ^5 of the form 



WW 1 



r.3 



r 



sin^(6') (89) 



If we take a general linear combination of these three in the form 



sm^{9)\Ar^ +Br + -\ (90) 



r 



we obtain a valid stream function. The particular choice 

vl.,.-sm^(^)|.^- — + -| (91) 

gives the well-known stream function for very viscous flow around a sphere of radius a and uniform flow at 
at rate U at inflnity. Having non-dimensionalized we would, for example, scale so that a — 1. In general 
we have a contour integral technique for solving the PDE given by equation (4.7). We should note of course 
that this discussion is conflned to showing only that the contour integrals contain the well known historical 
solution for a sphere due to Stokes. The problems associated with this solution are also well known - see for 
example Chapter 8 van Dyke (1964) for a discussion of the Whitehead paradox and its resolution by matched 
asymptotic expansions. 
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5 Small but non- vanishing Reynolds number 

A natural question to ask is to wonder how much of the above is dependent on the purely linear structure 
that arises in the biharmonic limit? We cannot yet answer this question for a general Reynolds number 
in three dimensions, but we can observe that something very interesting happens when we (a) go back to 
two dimensions and (b) consider the case of a small but non-zero Reynolds number. Let us go back to the 
non-dimensional form of Eq. (2.13). This is 

In terms of the complex variable w — x + iy, we can write this in the form 

. ay _ R fdjj d^tp _ djj d^jj \ 
dw^dw^ 2 ydw dwdiju'^ dw diJbduP' J 

Rather than pursuing the approach of Legendre (1949) and Ranger (1991, 1994) we can consider instead the 
interesting physical case of small but non- vanishing Reynolds number. Let us assume that the solution for ip 
may be written as 

V; = Vo + i?V'i+0(i?2) (94) 

and that 

^o = ^{wfo{w)+go{w)] (95) 

This is a very strong assumption, and it is well known that this assumption of a power series dependence 
on the Reynolds number may fail. There may not indeed be a sensible form for V'o over a simple domain of 
interest. The reader is again referred to Chapter 8 of van Dyke (1964) for a discussion of the Stokes paradox 
for a cylinder in a uniform flow. Our purpose here is to illustrate that the low Reynolds number perturbation 
equation may indeed be integrated using holomorphic methods. The result may be of use in refining the 
results for a certain sub-class of problems where there is both a meaningful ^o o-nd the inertia terms in the 
Navier-Stokes equations (i.e. the non-linear terms) arising from -00 remain small over the entire domain of 
interest. Under these strong assumptions we can proceed. The equation for -01 is, under these assumptions, 

i-^J^ = \ilM[^fl,{w)+%i^)+l^i^-h^ (96) 

This may be solved almost explicitly as follows. We let F{w),G{w),H{w) be holomorphic functions with the 
properties 

F'{w) = foiw) , G'{w) - goiw) , H"iw) - foHfoH (97) 

Then a particular solution to Eq. (5.5) is given by 

1 / u;2 \ 

V>ip = - f {wF'{w) - 2F{w))F{w) + F'{w)G{w) + —H(w) j (98) 

and a complementary function exists in the obvious form 

V'iCF = 3ft;{w/i(w)+5i(w)} (99) 

where /i and gi are arbitrary holomorphic functions. So we can see that apart from the practical issue on 
constructing the integrals in Eq. (5.6), the first perturbation can be constructed by separate integration of 
the w and w components. In fact, we have shown that the perturbativc non-linear problem may be solved in 
terms of free holomorphic functions F,G, fi, gi and the solution, apart from the construction of H, is given 
explicitly in terms of this holomorphic information. This observation gives some hope that a corresponding 
three-dimensional structure might exist. 
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6 Suminary 

We have presented the theory of 3D twistors in such a way as to ahow the understanding of biharmonic flow in 
three dimensions. The representations developed are a natural generalization of the stream function method 
that is very familiar to fluid dynamicists. Hopefully this work will stimulate the use of twistor methods in 
dealing with 3D problems in applied mathematics and fluid dynamics in particular. The ultimate goal is a 
better understanding of the Navier-Stokes equations through the geometry of complex variable techniques. 
Ideally we would be able to draw together the threads of work given by Legendre (1949), Ranger (1991,1994), 
and more recently, Roulstone et al (2005). Another major issue is the imposition of boundary conditions 
directly in twistor space. But it is quite clear that we can generate elegant contour integral solutions to 
problems in fluid dynamics, beyond those already well known for two dimensions and potential flow. We have 
also seen that the case of small but non-zero Reynolds number is open to complex variable treatment in the 
case of two dimensions, though at present we need to make some strong assumptions on the nature of the flow. 
Finally the generality of the proposed representations needs to be established. So far we have shown how to 
construct an interesting class of solutions, but further work is needed on just 'how large' this class is. 
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